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SUMMARY/OVERVIEW: 

Studies  were  conducted  to  finalize  a  fundamental  and  predictive  reaction  model  for  the 
combustion  of  hydrogen  and  carbon  monoxide  through  a  multi-parameter  optimization.  These 
studies  showed  that  reliable  data  of  hydrogen  and  carbon  monoxide  oxidation  at  high 
temperatures  can  be  reconciled  by  a  single  kinetic  model.  An  advanced  approach  to  solve  the 
master  equation  of  collision  energy  transfer  was  developed  to  predict  the  rate  constants  of 
combustion  reactions  of  arbitrary  complexities.  Studies  were  also  conducted  on  a  class  of  unique 
combustion  reactions  that  involve  spin  state  crossing — one  of  the  last  few  unresolved  theoretical 
problems  in  combustion  reaction  kinetics. 

TECHNICAL  DISCUSSION: 

1 .  Optimization  of  H>/CO  Combustion  Chemistry 

In  this  collaborative  project  (with  Professor  Fokion  Egolfopoulos  of  USC),  we  propose  a  H2-CO 
kinetic  model  which  incorporates  the  recent  thermodynamic,  kinetic,  and  species  transport 
updates  relevant  to  high-temperature  H2  and  CO  oxidation  [1].  Attention  has  been  placed  on 
obtaining  a  comprehensive  and  kinetically  accurate  model  able  to  predict  a  wide  variety  of  H2- 
CO  combustion  data.  The  model  was  subject  to  systematic  optimization  and  validation  tests 
against  reliable  H2-CO  combustion  data,  from  global  combustion  properties  (shock-tube  ignition 
delays,  laminar  flame  speeds,  and  extinction  strain  rates)  to  detailed  species  profiles  during  H2 
and  CO  oxidation  in  flow-reactor  and  in  laminar  premixed  flames.  Although  some  of  the 
features  of  model  optimization  have  been  reported  in  the  last  contractors  meeting,  the  model 
presented  herein  represents  our  final  effort  in  this  area.  Validation  tests  now  include  the 
extinction  stretch  rates  of  ultra-lean  H2-air  flames,  not  considered  previously  [2].  The  optimized 
model  will  be  used  as  a  kinetic  foundation  for  future  studies  of  aliphatics  combustion 
mechanisms. 

The  unoptimized  (trial)  reaction  model  consists  of  14  species  and  30  reactions.  The  trial 
model  was  based  on  a  careful  review  of  recent  kinetics  literature,  considering  both  direct  data 
and  compilations.  The  diffusion  coefficients  of  key  pairs  have  been  updated,  as  reported  in  the 
last  year’s  annual  report.  These  updates  are  reported  in  [2].  A  comprehensive  review  was 
conducted  for  a  large  number  of  H2-CO  combustion  data.  Thirty-six  experiments  were  chosen  as 

20050218  077 


optimization  targets.  The  targets  can  be  classified  into  four  categories:  (1)  laminar  flame  speeds 
of  H2-air,  Hh-CVHe,  and  H2-CO-air  mixtures,  (2)  the  peak  mole  fractions  of  H  and  O  in  a  low- 
pressure  burner-stabilized  Hh-CVAr  flame,  (3)  the  consumption  rates  of  H2  and  CO  during  the 
reaction  of  H2-O2-N2  and  CO-O2-H2O-N2  mixtures  in  a  turbulent  flow  reactor,  and  (4)  ignition 
delay  times  of  H2-02-Ar  and  H2-C0-02-Ar  mixtures  behind  reflected  shock  waves.  Active 
parameters  were  identified  through  sensitivity  analysis.  The  entire  set  of  active  parameters 
consists  of  28  ^-factors  and  third-body  efficiency  factors.  Minimization  was  carried  out  on  the 
objective  function  L 2  =  2,[(  Tfoxpt- Tfrcatc)/ 0/]2  subject  to  the  constraints  that  the  active  parameters 
were  bound  by  the  ranges  of  their  respective  uncertainties.  Each  target  was  individually  weighted 
by  their  uncertainty. 

Selected  validation  results  of  the  optimized  model  are  briefly  discussed  here.  Fig.  1  presents 
experimental  and  computed  laminar  flame  speeds  of  Eh-air  and  air-equivalent  mixtures  where  N2 
was  replaced  by  Ar  or  He.  It  is  seen  that  the  flame  speeds  predicted  with  the  optimized  model  are 
in  an  excellent  agreement  with  the  experimental  data.  Figures  2  shows  that  the  optimized  model 
reproduces  experimental  H2-CO-air  laminar  flame  speeds.  Fig.  3  shows  a  plot  of  experimental 
and  computed  ignition  delay  times  for  H2-C>2-Ar  mixtures  behind  reflected  shock  waves.  Here 
the  experimental  shock-tube  ignition  delay  data  were  fitted  into  r(ps)=[H2]"0154[O2]0  693[Ar]004 
[6.77x1  O'8  7°  252  e9234/r]  for  non  “run-away”  data,  i.e.,  those  fall  on  the  nearly  linear  portion  of  the 
curves  of  Fig.  3,  where  []  denotes  concentration  in  mol/cm3.  Figure  4  depicts  species  profiles  for 
4  selected  H2  oxidation  experiments  in  a  turbulent  flow  reactor.  Similar  results  for  moist  CO 
oxidation  are  presented  in  Fig.  5,  over  an  extended  pressure  range.  The  trial  and  optimized 
models  were  also  tested  against  extinction  strain  rates  of  ultra-lean  H2-air  mixtures  [23. 

The  optimization  procedure  also  allows  us  to  probe  the  residual  kinetic  uncertainties.  For 
example,  we  found  that  the  rate  coefficient  of  H+02=0H+0  always  stayed  within  5%  of  the  trial 
value  for  all  optimization  runs  made,  including  the  use  of  a  smaller  number  of  targets  and/or  a 
reduced  dimensionality  of  active  parameter  space. 

This  work  is  reported  in  Refs.  1  and  2. 

2.  A  Monte  Carlo  Master  Equation  Solver 

The  temperature  and  pressure  dependence  of  the  rate  constants  of  chemically  activated 
reactions  is  essential  in  combustion  modeling.  Many  gas  phase  reactions  proceed  through  the 
formation  of  ro-vibrationally  excited  adducts.  The  fate  of  the  excited  species  is  determined  by 
several  competing  processes.  Energy  is  transferred  to  and  from  the  molecule  as  a  consequence  of 
inelastic  collisions  with  the  surrounding  bath  gas  molecules.  Between  collisions,  the  molecule 
can  undergo  internal  rearrangement  or  isomerization,  decompose  to  form  new  products  or  simply 
dissociate  back  to  the  reactants.  Much  progress  has  been  made  in  our  understanding  of 
collisional  energy  transfer  and  its  influence  on  the  temperature  and  pressure  dependence  of 
unimolecular  and  chemically  activated  reactions.  Troe  introduced  a  collision  efficiency  factor 
J3C. ,  such  that  the  rate  constant  at  low  pressure  was  related  to  the  strong  collision  rate  constant 
through  k0  =  pji* .  The  modified  strong  collision  (MSC)  approach  thus  allowed  for  weak 
collisions  with  the  efficiency  of  energy  transfer  dependent  on  the  nature  of  the  bath  gas.  Though 
the  modified  strong  collision  model  was  intended  for  use  only  for  single-channel  unimolecular 
reactions,  over  the  years  this  model  has  been  widely  used  for  multi-channel  chemically  activated 
reactions.  The  fundamental  validity  of  this  approach  is  certainly  questionable.  For  this  reason  a 
Monte  Carlo  code  was  developed  that  can  solve  the  master  equation  of  collision  energy  transfer 
for  unimolecular  and  chemically  activated  reactions  of  an  arbitrary  complexity.  The  code  was 
fully  validated  and  optimized.  Comparisons  with  the  MSC  model  were  made.  These 
comparisons  show  that  the  MSC  model  lead  to  large  errors  in  rate  constant  calculations. 


A  second  motivation  for  the  development  of  the  Monte  Carlo  code  arises  from  our  interest  to 
tackle  a  fundamental  kinetic  problem  not  tackled  previously,  namely  chemically  activated 
reactions  that  involve  spin-state  crossing.  A  methodology  is  currently  being  developed  and  tested 
for  these  reactions,  including  those  to  be  discussed  in  section  3. 

This  work  is  reported  in  Ref.  3. 

3.  Reactions  on  C7H4O  and  CfiFhO  Potential  Energy  Surfaces. 

There  are  several  motivations  for  this  work.  The  reactions  on  the  C2H4O  and  CgHeO 
potential  energy  surface  involve  multiple  triplet  and  singlet  spin-state  crossings.  These  reactions 
are  ideal  for  use  as  test  cases  of  our  Monte  Carlo  approach.  The  reactions  on  the  CeHgO 
potential  energy  surface  are  critical  to  the  fundamental  kinetics  of  aromatics  oxidation.  The 
reaction  model  of  aromatics  is  currently  under  development  by  our  research  group.  The  C2H4O 
reactions  are  important  to  an  predictive  ethylene  oxide  reaction  model  also  under  development. 
The  letter  effort  was  initiated  following  the  request  by  the  USAF  Research  Laboratory. 

Quantum  chemistry  calculations  at  G3  and  CASSCF  levels  of  theory  were  carried  out.  The 
potential  energy  surfaces  for  both  C2H4O  and  C6H60  are  now  complete,  including  all  feasible 
spin-state  crossing.  Unique  kinetic  features  resulting  from  spin-state  crossing  will  be  examined 
using  the  developed  Monte  Carlo  code. 

Parts  of  this  work  are  reported  in  Ref.  4. 

TECHNOLOGY  TRANSFER 


The  following  technology  transitions  and  transfers  were  made: 


Performers: 

Customer: 

Phone: 

Result: 

Application: 


Dr.  Hai  Wang,  University  of  Delaware,  (302)  831-4464 

Dr.  Fokion  Egolfopoulos,  University  of  Southern  California,  (213)  740-0480 

Dr.  Bala  Varatharajan,  GE  Global  Research  Center,  One  Research  Circle, 

Building  K-l,  Room  ES243,  Niskayuna  NY  12309. 

(518)387-4133 

Reaction  model  for  hydrogen  and  propane  combustion 
NOx  control  in  gas  turbines 


Performer: 

Customer: 

Phone: 

Result: 

Application: 


Dr.  Hai  Wang,  University  of  Delaware,  (302)  831-4464 

Dr.  Biswadip  Shome,  GE  Gas  Turbines  (Greenville)  LLC,  300  Garlington  Road, 
P.O.  Box  648,  Mail  Sorter  GTTL  1256,  Greenville,  SC  29602-0648. 

(864)  254-4623 

Reaction  model  for  hydrogen  and  carbon  monoxide  combustion 
Flame  holding  in  gas  turbines 


EXTRAMURAL-INTRAMURAL  INTERACTIONS 


We  worked  closely  with  Dr.  Dean  Eklund  and  Dr.  Thomas  Jackson  until  their  departure  from  the 
Wright-Patterson  AFRL.  Currently,  our  AFRL  contact  person  is  Dr.  J.  T  Edwards.  A  visit  to  the 
Wright-Patterson  AFRL  has  been  scheduled  for  September  15,  2004.  During  the  visit,  the  PI 
will  meet  and  discuss  potential  research  collaborations  with  various  lab  research  personnel,  and 
give  a  talk  entitled  “Incipient  Soot  Formation  and  Growth.” 
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Figure  1.  Experimental  (symbols)  [Tse,  S.  D.,  Zhu,  D.  L.,  and  Law, 
C.  K,  Proc.  Combust.  Inst.  28,  1793-1800  (2000);  Taylor,  S.  C., 
"Burning  Velocity  and  the  Influence  of  Flame  Stretch,"  Ph.D.  thesis. 
University  of  Leeds,  1991;  Aung,  K.  T.,  Hassan,  M.  I.,  and  Faeth,  G. 
M.,  Combust.  Flame  109, 1-24,  (1997);  Egolfopoulos,  F.  N.  and  Law, 
C.  K.,  Proc.  Combust.  Inst.  23,  333-340  (1990);  Kwon,  O.  C.  and 
Faeth,  G.  M.  Combust.  Flame  124  590-610  (2001)]  and  computed 
(line)  laminar  flame  speed  of  H2-air  and  H2-Ar-  and  He-substituted  air 
mixtures  at  the  atmospheric  pressure. 


Figure  2.  Experimental  (symbols)  {Vagelopoulos,  C.  M.  and 
Egolfopoulos,  F.  N.  Proc.  Combust.  Inst.  25,  1317-1323  (1994)]  and 
computed  (line)  laminar  flame  speed  of  H2-CO-air  mixtures  at  the 
atmospheric  pressure. 


delay  times  of  H2-02-Ar  mixtures  behind  reflected  shock  waves. 
Symbols:  (a)  6.67%H2+3.33%02,  />5=1 .3  5-2.90  atm;  (b) 

5%H2+5%02,jp5=l  .35-2.90  atm  (onset  of  pressure  rise,  [R.  K.  Cheng, 
A.  K.  Oppenheim,  Combust  Flame  58,  (1984)  125-139]),  (c) 
0.5%H2+0.25%02,  p$=33  atm,  (d)  2%H2+i%Q2,  py=33  atm,  (e) 
0.5%H2+0.25%02,  P5=57atm,  (f)  0.33%H2+0.17%02,  ps=6 4  atm,  (g) 
0.1%H2+0.05%02,  p5= 64  atm,  (h)  0.5%H2+0.25%02,  pf=%l  atm 
(maximum  OH  absorption  rate,  [E.  L.  Peterson,  D.  F.  Davidson,  M. 
Rohrig,  R.  K.  Hanson,  20th  lnt.  Symp.  on  Shock  Waves,  1996,  pp.  941- 
946]),  (i)  8%H2+2%02,  p$~5  atm  (maximum  OH  emission)  [G.  B. 
Skinner,  G.  H.  Ringrose,  J.  Chem.  Phys.  42  (1965)  2190-2192],  (j) 
four  H2+02  mixtures  [A.  Cohen,  J.  Larsen,  Report  BRL  1386,  1967]. 
Lines:  (1)  0.5%H2+0.25%02,  p5=Sl  atm  (maximum  [OH]  rate),  (2) 
2%H2+1%02,  p$-33  atm  (maximum  [OH]  rate),  (3)  8%H2+2%02, 
p5-S  atm  (maximum  {OH]),  (4)  5%H2+5%02,  p$=2  atm  (maximum 
pressure  gradient). 


Figure  4.  Experimental  (symbols)  [Mueller,  M.  A.,  Kim,  T.  J., 
Yetter,,  R.  A.,  and  Dryer,  F.  L.  lnt.  J.  Chem.  Kinet.  31,  113-125 
(1999)]  and  computed  (lines)  species  mole  fraction  profiles  during 
hydrogen  oxidation  in  a  flow  reactor.  Solid  lines:  optimized  model; 
dashed  lines:  trial  model. 


Time  (msec) 

Figure  5.  Experimental  (symbols)  [Kim,  T.  J.,  Yetter,  R.  A.,  Dryer,  F. 
Proc.  Combust.  Inst.  25,  759-766  (1994)]  and  computed  (lines) 
(CO]/[CO]o  profiles  during  moist  CO  oxidation  in  a  flow  reactor. 
Cases  (a):  L014%CCH-0.517%02  +0.65%H2O  in  N2,  p=  1  atm, 
r0=1038  K,  (b)1.01%CO  +0.496%02+0.65%H20  in  N2,  p=2.44  atm, 
1038  K,  (c)  0.988%CO+0.494%02+0.65%H20  in  N2,p=3.46  atm, 
r0=1038  K,  (d)  0.984%C00.497%02+  0.65%H2O  in  N2, p=6.5  atm, 
T0=  1040  K,  (e)  0.994%CO  +I.47%O2+0.65%H2O  in  N2,  9.6  atm, 

7o=1039  K. 


